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Abstract: Towards modelling the charge asymmetry observed in heavy ion collisions, we 
present here analytic solutions of relativistic hydrodynamics containing parity violating 



and anomalous terms at the first order in the hydrodynamic approximation. These terms 
can induce chiral magnetic and chiral vortical effect leading to the generation of the charge 
asymmetry. We also consider sphaleron solutions with non trivial winding number to model 
the phenomenon. We calculate the net chiral charge difference produced in our solutions. 
We anticipate their relevance also in the context of baryogenesis in early universe, neutron 
star and some condensed matter situations. 
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1 Introduction and motivation 

For the early phase during the formation of quark gluon plasma in heavy ion collisions, 
relativistic hydrodynamics offers a good model to describe the dynamics. Relativistic 
hydrodynamics extends the regime of applicability of hydrodynamics to the situations of 
fast moving fluids such as plasma. Such cases are also encountered in nuclear physics and 
astrophysics apart from the heavy ion collision experiments. The equations governing the 
dynamics of the relativistic fluid are conservation equations for the stress energy tensor 
and the conserved current (s). These are relativistic analogs of the continuity equation and 
Navier Stokes equations. 

Recently, it was shown that the relativistic hydrodynamics can contain terms which 
do not have any analogs in non-relativistic cases. The conserved current can contain parity 
violating terms proportional to vorticity and/or magnetic field. These terms were first 
noticed during the investigation of equations governing small perturbations at the boundary 
of the charged black branes and drawing their parallels with Navier-Stokes equations. [1, 
2, 3] The theory of gravity in such cases is related to the strong t' Hooft coupling limit of 
a dual large N field theory and only the long range modes survive in the hydrodynamic 
approximation, i.e. modes surviving in long temporal and large wavelength limit. These 
extra terms were also later understood in terms of quantum triangle anomalies. [4] The 
anomalies are usually calculated in the perturbative limit of the quantum theory where 
the coupling is small. The topological nature of anomalies has a role to play to account 
for their occurrence in both the approaches to the quantum theory. We will discuss here 
two recently discovered terms in the expression of the current which are also related to 
anomaly of the previously conserved current. Both of them are parity odd terms. First one 
is proportional to vorticity and it leads to the phenomenon of chiral vortical effect. The 
second term is proportional to the magnetic field and results in chiral magnetic effect. [5]. 
For non abelian hydrodynamics particularly suited for quark gluon plasma, an effective 
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action for fluid with anomalies was constructed and chiral magnetic effect was obtained 
from it in [6]. 

In the case of quark gluon plasma generated during collision of heavy ions, the anoma- 
lous current of interest can be axial current leading to net chirality difference along the 
direction of background magnetic field. The net chirality difference can induce an electric 
field and a flow of charged carriers along the direction of magnetic field. This is known 
as chiral magnetic effect. [5] The rate of chiral charge difference depends on the strength 
of the anomaly. For a fluid with certain configurations for its velocity field (it'*), similar 
buildup of the electric field along the vorticity vector (u^ = ^e^ up<J u u d p u a ) happens during 
the process of chiral vortical effect. Such processes are not only restricted to quark gluon 
plasma, but are also important in many other contexts (to be discussed below). In all 
these cases, relativistic hydrodynamics offers itself as a good candidate model to describe 
the dynamical process. 

In quark gluon plasma, the gluon field can locally have some topologically non triv- 
ial configurations. These configurations are classified by their winding numbers and are 
separated from each other by a potential barrier of the order of Aqcd- The system can 
tunnel through it due to an instanton. But such processes are exponentially suppressed 
at high densities and at weak couplings. [7] The system can also roll over the barrier and 
such transitions are called sphalerons which manifest themselves during the presence of 
anomalies. These metastable configurations are quantum fluctuations possibly of the order 
of few fms and lead to a difference in left and right handed quarks. The discrete symmetries 
such as parity and CP are locally violated in these spatially localized bubbles. A sphaleron 
in non-abelian gauge fields in QCD can offer an explanation to the strong CP problem. 
They can also be present in quark-gluon plasma in a high magnetic field. [8, 9] Thus, they 
can lead to an observable charge separation by generating the net chiral difference and 
subsequent interaction of chiral particles with background magnetic field leading to chiral 
magnetic effect. 

The sphaleron offers a mechanism to explain many different effects in different settings. 
For example, sphaleron processes can be partly responsible for the baryon asymmetry in 
our universe. Before electroweak symmetry breaking, a non trivial topological configura- 
tion of hypercharge electromagnetic field could have occurred in the early universe plasma. 
Such a sphaleron configuration offers a model to generate the requisite amount of baryon 
asymmetry needed for the subsequent process of nucleosynthesis. [10] These processes vio- 
late baryon number conservation as well as local P and CP symmetries. The rate of baryon 
asymmetry production was estimated in [11] and it was found that such explanations are 
favourable if the electroweak symmetry breaking in the early universe was of first order. 
Recently, early universe baryogenesis due to anomaly in lepton number current was also 
considered. [12] 

The anomalies in 4 dimensional field theories can also be understood in terms of their 
dual gravity solutions (if they exist) in the context of AdS/CFT correspondence. [13, 14, 15, 
16, 17] The well understood N=4 Super Yang Mills theory can have anomalous R charge 
currents. It maps to a nontrivial 5 dimensional bulk Chern-Simons term in the gravitational 
Lagrangian. The magnetic field acting as a source in the chiral magnetic effect amounts 
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to putting appropriate boundary conditions on the bulk electromagnetic tensor F^ u in the 
dual gravity solution. The anomalous current will correspond to another abelian gauge field 
in the bulk. Holographic computations were also performed to determine chiral magnetic 
effect for the case of anisotropic fluid and its dependence on elliptic flow coefficient t>2, a 
quantity which parameterizes the event charge anisotropy in heavy ion collisions. [18]. 

Superfluids can also display a similar phenomena called chiral electric effect. [19, 20] 
Here, an anomalous current is generated due to a topologically non-trivial configuration 
of electric field. Moreover, chiral vortical effect in a pionic superfluid leads to a flow of 
fermionic zero modes along the direction of vorticity.[21] For chiral magnetic effect in the 
same medium, there will be strings carrying magnetic flux and the magnetic field plays the 
role of the vorticity vector. Chiral magnetic effect can also be present in metal crystals 
having a non trivial Berry phase configuration in the presence of electromagnetic field. [22] 
If there are k quanta of Berry curvature flux associated with any given Fermi surface, then 
the fermionic number current will be anomalous with the anomaly proportional to kE.B. 
This triangle anomaly will also give its contribution to the density-density correlator. If 
there are 2 Fermi surfaces with unequal chemical potentials in the presence of magnetic 
field, the chiral magnetic effect will manifest itself as a flow of fermionic current between 
the Fermi surfaces with a strength proportional to the magnetic field and the difference in 
chemical potentials. 

In this paper, we will explore these parity violating effects in the case of quark gluon 
plasma using the hydrodynamic approach. We will partly develop and demonstrate two 
methods of constructing hydrodynamic solutions containing these parity violating and 
anomalous terms. We will keep the dissipative coefficients vanishing throughout, so our 
fluid can be interpreted as a perfect fluid in the presence of parity violating and anomalous 
coefficients. In quark gluon plasma, the viscosity is actually very small. [23] In section 2, 
we write down the equations for the relativistic hydrodynamics at the first order explain- 
ing our conventions. Some solutions witnessing chiral vortical effect can be constructed 
by first finding a relativistic generalization of some known non-relativistic solutions and 
then modifying them to admit the non-trivial vorticity terms. This is demonstrated in 
section 3 for the case of a famous non relativistic solution known as Taylor-Green vortex. 
We will keep the electromagnetic fields vanishing here and calculate the net axial charge 
difference generated which amounts to chiral vortical effect. The second method is to use 
Hopf mapping to construct a topological solution with winding number one using the ve- 
locity and electromagnetic fields. A non relativistic magnetohydrodynamic solution based 
on such mapping was constructed earlier in [24]. However, we find that Hopf mapping 
can be used to generate a larger set of solutions of relativistic hydrodynamics, potentially 
setting a stage to explore in detail many dynamic processes. This will be the content of 
section 4. We will also find some simple solutions using this method in this section. In 
section 5, we use the same method to generate a sphaleron solution in relativistic hydro- 
dynamics with a topologically non trivial configuration of the background electromagnetic 
field. This solution has all the parity violating and anomalous coefficients non trivial at 
the first order. Even though, topological configurations for non-abelian fields seem to be 
more interesting, we will restrict ourselves to U(l) fields only in this paper for simplicity. 
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In the case of quark gluon plasma, it can be thought of as a restricted abelian version of 
chromo-electromagnetic fields or a background U(l) field produced by the highly energetic 
colliding charged particles. Also, all the quantities denoting space-time dimensions are kept 
dimensionless throughout in this paper. The point of view is to assume some natural length 
scale present in the theory and then dimensionless position variables denote multiples of 
it. The natural length scale in the case of quark gluon plasma can be inverse of either 
dynamically generated QCD scale, Aqcd or center of mass energy, s cm or temperature, T . 
We consider our solutions as representative solutions displaying chiral effects and calculate 
chiral charge difference contributed by them. Furthermore, we relax the boundary con- 
ditions to simplify the construction of the analytic solutions. In many cases of potential 
applications as mentioned above, boundary conditions may not play any significant role 
in determining the underlying processes. We also didn't impose the equation of state and 
it is determined implicitly by the additional assumptions that we make to simplify our 
equations. We consider physically interesting cases to be those for which the pressure and 
energy density are both positive everywhere. 

2 Equations for anomalous hydrodynamics 

The equations for relativistic hydrodynamics are given as conservation equations for the 
stress-energy tensor, T^ u and the conserved current(s). See eqs. (2.1). We will consider 
the case of one conserved current, j M which can represent particle flux like baryon current 
or axial current and is likely to get anomalous contributions. We also have background 
electromagnetic fields with the electric and magnetic fields defined as = F^ u u v and 
= he^ va ^u v F a p, where F^ v is an antisymmetric field strength tensor. In the fluid 
rest frame, their spatial components indeed represent electric and magnetic fields. The 
conservation equations are supplemented with constituent equations which express stress 
energy tensor(T' 1I/ ) and current^) in terms of pressure(P) , enthalpy density(Zi), tem- 
perature (T), particle number density (n), velocities (u M ) and their derivatives. Fluid is 
likely to have non trivial vorticity defined as = i2e^ upa u v d p u a giving rise to a parity 
violating term in the current conservation equation. In the presence of electromagnetic 
fields, one can have a parity violating term proportional to magnetic field as well as an 
anomalous term in the current equation. The presence of such terms was also understood 
in terms of constraints on near equilibrium partition function. [25]. They were also calcu- 
lated using Kubo formula and do receive corrections proportional to gravitational anomaly 
coefficient. [16, 26]. We write below equations for relativistic hydrodynamics containing 
these terms. 



= F vX j x 

d^f = —e^F^Fpv = CE^B^ 

= hu»u v + PgT - V P> ta P v f > (d a Uf ) + d p u a ) -((- ^v) P^dxu 
f = nu » - aTP^du (|) + aE" + ^ + ( B B^ 
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These are supplemented with the relativistic constraint on velocity, u^u u = —1. Here, 
H is the chemical potential and g^ v is the metric, which we take to be Lorentzian with 
signature (—,+,+,+). The notation, P nv = g nv + u n u u denotes the projection tensor 
and projects any tensor perpendicular to the velocity field. It is a symmetric tensor and 
satisfies relation u n P^ u = 0. The density n may represent axial charge density in the case 
of quark gluon plasma or baryon/lepton charge density in the case of baryogenesis in the 
early universe. The various dissipative coefficients are bulk viscosity (£), shear viscosity 
[rj) and conductivity (a) . We will call the coefficients of parity violating terms £ and £b as 
chiral vortical conductivity and chiral magnetic conductivity, respectively. The coefficient 
C denotes the strength of the anomaly. It is related to chiral conductivities as 

U-c{,-\^) (2,) 

The electromagnetic fields also satisfy Maxwell equations i.e. field strength conservation 
equation and Bianchi identity. 

V = Jem, 

The quantity j EM represents the background electromagnetic four-current i.e j EM repre- 
sents the charge density and its spatial part, j l EM represent the electric current, where i = 
1, 2, 3. The convention for Levi-Civita symbol used in this manuscript is eoi23 = e" 1 12 = !• 



3 Solution by uplifting 

In this section, we will try to estimate the chiral effects due to a fluctuation plausible in 
quark gluon plasma. To do so, we will first construct an explicit relativistic hydrodynamic 
solution by uplifting a solution of non-relativistic hydrodynamics. Even if we expect that 
the observable quantities in quark gluon plasma may not strongly depend on the details of 
such fluctuations, construction of few explicit representative solutions will make possible 
any explicit computations of relevant quantities. Here, we choose a well known solution of 
non-relativistic hydrodynamics called Taylor-Green solution to estimate its chiral effects. 
The non-relativistic equations are given as 

dv 1 ^ dv 2 _ q 
dx dy 

dv 1 .dv 1 ndv 1 1 dP n ( d 2 v 1 d 2 v 1 
+ v 1 — + v 2 — = + v —r + 



dt dx dy p n dx \ dx 2 dy 2 

dv 2 l dv 2 1 ®Pn ( d 2 v 2 d 2 v 2 \ 

dt dx dy p n dy \ dx 2 dy 2 J 

Here, v i (i = 1,2), P n , e n , and p n denote velocity components, pressure, energy density 
and mass density of a non-relativistic fluid in 2+1 dimensions. They are continuity equation 
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and conservation of momentum flux for constant mass density p n . We then try to modify 
their solutions to also accommodate new coefficients. The first part can be done using the 
prescription given in [27, 28]. According to it, given a solution of non-relativistic equations 
in 2+1 dimensions, a solution of relativistic equations in 3 + 1 dimensions can be written 

as 



u 



u 



1 



Pr. 



2 £n Pn 



1 



+ v/vT 



u 

p 



u + v l 



Pn 



p = 2e n + p r , 



(3.2) 



Here, u + and u~ denote the velocity components along the null directions i.e. = 
-^(u° ±u z ) and u l , the same along other two spatial directions. The coordinates along the 
three spatial directions are denoted as x, y, and z. Symbols P and p denote relativistic 
pressure and density of the fluid. 

The Taylor Green vortex solution is given as in [29], 



v 

F(t) 

Pn 



F{i) sin x cosy, 
—F(t) cos x sin y, 



-2ut 



^F(t)(cos 2x + cos 2y). 



(3.3) 



We consider the simple case of zero viscosity. We put v = here and then try to get 
the relativistic version of it. This solution has a relativistic analog in 3+1 dimensions as 
discussed in eqs. (3.2). 

en 1 1-1/2 
2— H — (cos 2x + cos 2y) 

Pn 2 



u 



1 + 



2er 

Pn 



2 sin xsin y 



u 



u sin x cosy, 
— u + cos x sin y, 



P 
h 



■(cos 2x + cos 2y), 



(u+) 2 ' 

The above quantities satisfy the relativistic equation for perfect fluid i.e. 



0. 



(3.4) 



(3.5) 



To find the modified solution of relativistic equations with non-trivial new coefficients, we 
first simplify the relativistic hydrodynamic equations. We choose the first order dissipative 
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coefficients £ = a = r/ = and write F^ v in terms of electric and magnetic fields as 

Ffiv = 2u {ll E v] - e^ pa uPB a (3.6) 

The conservation equations are then 

= nE v + u^bE^ + ZE p ljP - tu pZ P] - iz v 

d„f = CE^ (3.7) 

where 

z» = B a dK a ^ (3.8) 

We then simplify the equations by taking the background electromagnetic fields to be 
vanishing, i.e. Ea = B p = 0. We further break the equation with stress tensor into its 
trace and traceless part by taking the trace with the velocity vector, u^. After a little 
algebraic manipulation, we get the following equations, 

n3 M (W) - hf(u.d)u p = (j.d)P, 

d.j = o, 

hu^d p u u + u u u^d p P + d u P = 0. (3.9) 

We find that the relativistic analog constructed in eqs. (3.4) helps a lot in dealing with these 
equations. To accommodate the non trivial vorticity term with chiral vortical conductivity 
£, it suffices to choose a suitable expression for the number density n. We present our 
solution in terms of Lorentz factor x, which in this case turns out to be 

X -2 = — + -(cos2x + cos2y) (3.10) 

Pn 2 

This variable is equivalent to u + in eqs. (3.4). The solution can then be written as 

h 



Pn 
2 ' 



XV 11 , 



X 

u» = 

v + = 1, 
- 1 . e n -2 -2 

v = — I sin x sin y, 

2 Pn 

v 1 = sin x cosy, 

v 2 = — cos re sin y, 

P = P + ^(cos 2x + cos 2y), 

n = — xsinxsiny. (3-11) 
3 

We consider Pq to be a constant positive quantity. The chiral vortical conductivity appears 
explicitly in the last expression for the number density n. The only connection between 
the two conservation equations is the velocity. The current is 

— s 2 • ■ h , 2e n . 2 -2 

j = — x sin x sin y 1 H sin x — sin y 

3 V Pn 

3 + =f=j y = 0. (3.12) 
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Figure 1. Plot of maximum axial charge separation (in units of £L) as a function of ratio of energy 
and number density, e n / p n generated for the case of relativistic Taylor Green vortex solution. 



This solution is in a steady state. This is expected as we have dropped all dissipative terms. 
The non-trivial chiral vortical conductivity £, does not lead to dissipation. We calculate 
the zeroth component of vorticity to be 



X 

3^/2 



sin x sin y 



e n 3 2 2 
h cos x cos y 

Pn 2 



(3.13) 



The contribution to axial charge difference will be f d^x^u . We assume that the solution 
holds for some length L along the z direction. Along other two spatial directions, the 
contributions from different regions tend to cancel due to sinusoidal dependence. So, the 
net contribution depends on how far the solution extends. The maximum contribution will 
be for the situation in which the solution extends from to ir along both x and y directions. 
Then the axial charge difference generated will be 



a € L f n j r j sinxsiny(e n /p n - 3/2 + cos 2 x cos 2 y) 

Aha = — 7= / ax ay ; = 5 — . 

3V2 Jo Jo {Ztn/Pn - 1 + cos^ x + cos z y) 



(3.14) 



The axial charge difference created will induce an electric field resulting in chiral vortical 
effect. It is evaluated numerically and displayed as a function of e n /p n in fig. (1). 



4 Solution using Hopf fibration 

The background electromagnetic fields generated in quark gluon plasma are due to charges 
of the colliding ions. Due to the high energy present in the collisions, a non trivial con- 
figuration of the background field is also likely to occur. In this section, we try to find 
relativistic solutions with topologically non trivial background electromagnetic field. An 
interesting non-relativistic solution is given in [24], which has a non-trivial index defined 
as 

I = — / A(V x A)d 3 x. (4.1) 

The solution was obtained there using Hopf fibration map. Given such a map / : — > S 2 , 
one can pullback the volume form of 2-sphere to get a two form on 3-sphere. Any two form 
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on a 2-sphere can be written as an exact differential of a one form. The vector dual of such 
a one form can be projected on B? using stereographic projection. The value of the above 
index is one for the vector potential in E? generated by this procedure. The latter is given 
as [24] 

.1 _ (xz-y) 



2r 2 

2 _ (yz + x) 



2r 2 

■ 

A 



3 _ (2z 2 + 2 



4r 2 

where r = 1 + x 2 + y 2 + z 2 . (4.2) 

We will use this vector potential to generate solutions of relativistic hydrodynamic equa- 
tions which are more general than the simple relativistic generalization of the solution in 
[24]. We make the assumption that the profiles of the fluid and that of the background 
vector potential (A^) are proportional to the vector potential given above, i.e., 

A> 1 = {P,aA i ). (4.3) 

Here, /, v, f3 and a are functions of the radial coordinate r = 1 + x 2 + y 2 + z 2 and time t. 
The electro-magnetic tensor is calculated to be 



F ij = Qi^j _ Qj^i = € ijk 



i J • 



2a A vol 

a } A k + — a k 

r J * 

where =— (— y,x,l). (4.4) 

We denote the derivatives with respect to time and radial coordinate r by dot ' and prime 
', respectively. The Bianchi identity is automatically satisfied by this construction. We get 
the following expressions for the electric and magnetic fields. 



B 



o 



fa 
12r 3 ' 



B l = -| - vol + /?'/) At + ^(/37 - va'). (4.5) 

The vorticity here is 



24r 3 ' 

Ai (~2vf , , A ra. 



3 V r 



+ vf-v'f\+ r -^{v'f-vf). (4.6) 
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The equations that need to be satisfied are the conservation equation for stress energy 
tensor and the current equation along with the relativistic constraint on velocity. 



%TV = F» x j x , 

T» v = hu^u v + Pif u 

= -1. (4.7) 



Here, we have taken the conductivity and viscosities to be zero (er = rj = £ = 0). These 
equations for our case reduce to 



A 1 



(vhf) + —^(rhf 2 )' + nva 
2r z 



2P' - ^ + 2nv0 



mfa' £ , , 2 , 

— 2 ^ ■qi( av ) f ~ av f I + 3?b«/P 



+ 2^3 +2/3'^)(e/ + 2£ B a) 



(^ 2 ) + ^{Tvhf)'-P 



2r 2 
1 



48r 3 



(d + 8r/3'z)(-3rnf + 2fu/ + 4£ B t;a), 



1 + 



JL 

16r 2 



+ 



(4- 



They are coupled non-linear partial differential equations and in order to solve them, we 
make an ansatz that the factors proportional to A 1 , zA l , x 1 and z cancel out separately. 
The resulting equations can be written elegantly in terms of two new variables defined as 



M = -mf + -v(Zf + 2ZBa), 



N 



nv 



f 



24r 3 



(£/ + 2£ B a) 



— f M+^-{U + 2i B a). 
rf 3rf 



(4.9) 
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We thus obtain the following set of equations. 



(vhf) + Na = 0, 



a'M + ^a(fv' - vf) + ^-a{fP ~ vol) = (rhf 2 )', 

P ' + + 0N = °' 

12r 3 ' 
(rvhf)' = f3'M, 

M' + ^(jV - vf) + ^(f{3' - vol) + \ca0 = 0. (4.10) 

We next choose Coulomb gauge A = (3 = without loss of generality and for simplicity, 
we look for only steady state solutions. In other words, we assume that all the functions 
explicitly appearing in the set of equations above are time-independent. We denote rvhf 
by a constant k. This reduces the above set of equations to 

p' = — K(hfy 

(v)' + {rnf)a> = Y [{av) ' f " {av)f ' ] 

(m/)' = p{Zf + SB<x). (4.12) 

Here, v and / are related by the constraint v = 1 + jg^- The equations, even with the 
strident looking assumptions, admit a wide class of solutions. We consider two cases. 
Case I: C = 0. 

This is a perfect fluid case with no chiral vortical and chiral magnetic conductivity 
and is the simplest non trivial solution of this class. The eq. (4.12) tells that rnf = Nq (a 
constant). The other equation (4.11) can then be solved to obtain the following solution. 

Akrf 

a - 



Ak 

N 
rf 

^"-^ w/^ lf^ (4I3) 

Here, Po represents the pressure at the origin, x = y = z = 0. The function /(r) is left 
undetermined and can be any function which leads to well defined physical quantities i.e. 
pressure, energy density and enthalpy density. 
Case II: v = A/, A is a constant. 



- 11 - 



In this case, the spatial part of vorticity is proportional to spatial part of velocity or 
background vector potential. Velocity constraint and equation for pressure leads to 



/ 
P 



4r 



k 
48A 



1 - ^ + Pq. 



(4.14) 



Here, Pq is a positive constant denoting the pressure at the origin. We find the number 
density using equation (4.11) to be 

3 r 

The fact that rvhf = k can be used to evaluate the enthalpy density. 



n 



(4.15) 



h 



k 



Ar/ 



2 ' 



(4.16) 



Equation (4.12) leads to a consistency relation stating that = which implies a = 0. 
This is then a case of pure gauge and hence no electromagnetic fields. One gets the 4- vector 
and j 1 * for this case to be 



Co" 



24r 3 ' 3r 



0,0,0 



2£ 

3r : 



Near the origin x = y = z = i.e. r = 1, the solution simplifies to 



(4.17) 



n 



V16A 2 - 1 

3V16A 2 - l' 
fc(16A 2 - 1) 



A, 0,0, 



1 



16A 



P = Pn. 



(4.18) 



To see the asymptotic behavior of the solution at large distances r — > oo, we use spherical 
coordinates 9,4>, so that x = \/r — 1 sin 8 cos cfi, y = \Jr — 1 sin 9 sin <fi and z = yr — 1 cos 9. 
Then as r — > oo, 



u» - 




n - 


3r ' 


h - 


fcA 
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P - 


" jPo + 48A- 



1 

4A7 



1 

4A7 



(4.19) 
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From the pressure profile, we see that such configurations can happen around local de- 
pressions in pressure. As r increases, the pressure increases to a constant value. The 
fluid velocity, number density and enthalpy decreases as 1/r. The magnitude of velocity 
or speed is spherically symmetric at large r. This solution suffers a problem that energy 
density becomes negative for r larger than r ~ 48A 2 (1 + 48A;A/Po) _1 - So, solution can be 
relevant only for small r and should be dominated by some other solution at larger r. It 
can happen that some of the dissipative coefficients neglected here can prevent the solution 
from this problem. To calculate the contribution to the axial charge difference caused by 
the vorticity, we assume the solution to hold upto r = R. We thus obtain 



An A = jd*x&» = ~f 
2vr£ 



4vr£ f R y/T^ldr 



^(16A 2 r 2 - 1) 



3^A 
8vr£ 



sec" 1 VR (4.20) 

This axial charge difference will induce an electric current and hence, will lead to chiral 
vortical effect. 



5 Sphaleron Solution 

We revisit the eqs. (4.10) again and simplify them by choosing v = Xf and /3 = Xa. We 
also choose a = kf. These assumptions make the spatial part of vorticity, magnetic field, 
velocity and vector potential proportional to each other. We look here for steady state 
solutions only. Then, the set of equations needed to be solved is 

(rhf 2 )' = a'M, 

M = --XCa 2 + a , 

(hf 2 )' 

P' + LL^ + Xa'N = 0. (5.1) 

We solve top two equations above to obtain 

4r 



Vl6A 2 r 2 - 1 



f 



h = --XCk 6 - + a k— + 

9 r rj rj z 

n= 2 -K {C k 2 + t + 2ki B )-^ f . (5.2) 
Sr rj 



We simplify the equation for pressure to get 

2 



P - ^ ( { + 2Kb + gCfj ^ _ + ^. (5.3) 
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This results in the expression for pressure as 

P = P + f AA; U + + |cfc : 



tan 1 T~ 



As r — t- oo, the asymptotic behavior of various expressions above are 



k 2 ( 2 
h -> — [ao + bo - -CX 
r \ 9 



bp 
48r 3 ' 



7? 



2A 

— {Ck 2 +t + 2kt B )-avk 



P -> Pn 



1 



48H 



2A 



The associated background electric and magnetic fields are calculated to be 



\kf\ kf j kp_ 
32r 4Z '32r 4Z 8r 3 X 



kf 2 _4Xk f2Ai 



12r 3 ' 3r 

The vorticity and the anomalous current are 

f 2 2A 



0J r 



24r 3 ' 3r 



/A 



r \ 3 r r / 



3r v r \ 3 

The background electromagnetic current is found to be 



jg M = 2Afc{2(r-l)/" + 3/'} 



3Xkf 
128r 5 



{1 + 16A 2 r(r - 2)}, 



Jem — 2kA l 



4 (2f 



f ) - 2f + 2 



2/ 



/' + (r/7 + 



Hrf'Y xi 



kfA 1 
128r 6 

kX 2 fa i 
16r 3 



12 



2r 



+ 2(7 - 2r) 



4r 



+ 3(2-r) 



(1 + 32A 2 r 



fc/ 5 

2 9 r 6 



(1 + 32AV)4. 



(5.4) 



(5.5) 



(5.6) 



(5.7) 
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There are 6 unknown parameters, namely A, k, C, ao, bo, Pq along with 2 chiral conductivi- 
ties £ and £g ■ Parameter A is related to the ratio of fluid velocity to speed of light (usually 
denoted as /3) as A = (4r/3) _1 . Parameters k and C denote the strength of the electro- 
magnetic field and the anomaly, respectively. Constants ao, &o and Po are related to the 
axial charge density, enthalpy and pressure in the absence of non-trivial topological charge 
for this solution. We take these constants to be zero as we are interested in effects due to 
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Figure 2. Plot of enthalpy density (red curve) 
and 100 times the pressure (blue curve) vs r. 
The maximum of these physical variables oc- 
cur at the origin. 




0.05 0.10 0.15 0.20 

Figure 3. Energy Density is plotted against 
pressure. The energy density is more sensitive 
to changes in pressure at higher pressure, i.e. 
near the origin. 



topological charge. We can then simplify the chiral conductivities related by eq. (2.2) to 
obtain 

Cfi 2 



1 - 2fi/k 



C/x(l-/V(2fc)) 
63 = 1 - 2^/fc (5 - 9) 

Real values of / needs A > 1/4. We consider physically interesting cases as those for 
which P and e = h — P are always positive. All possible choices of parameters do not lead 
to physically sensible results. Below, we present one such choice, though other choices can 
also be feasible. In quark gluon plasma, fluid velocity can be very high up to fractions 
of speed of light. [30] We choose parameter A = 1 which implies that the maximum fluid 
velocity in the configuration is a quarter of the speed of light. Quantity C, related to 
anomaly in the corresponding field theory , should be some number of order unity. We 
choose it to be C = 1. Strength of the electromagnetic field and chemical potential vary 
a lot from one region in space-time to other. [5] We choose k related to electromagnetic 
field to be —1. Minus sign helps us to get physically sensible results for this configuration. 
We choose /i = 1/2 for the same reason. It is also a typical value of chemical potential 
for strong magnetic field. [5] For this choice, the radial dependence of enthalpy density and 
pressure as well as the implied equation of state are plotted in figures 2 and 3. 

The anomaly present in the right hand side of eq. (2.1) i.e. CE^B^, can be written 
as a total derivative of a Chern Simons current given as 

C 
7 

The Chern Simons charge of this sphaleron is 

y/r — ldr 
r(16A 2 r 2 - 1) 

(V4A + 1 + V4A - 1 - 4\/A) (5.11) 



K» = --e» vpa A v F pa . (5.10) 



N CS = J d 3 xK° = iCnk 2 J 
Cn 2 k 2 



VA 
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The rate of topological winding number changing transitions caused by the sphaleron so- 
lution will be proportional to the above charge, Ncs- I n the case of quark gluon plasma, 
it will contribute towards to the rate of production of chirality difference and the induced 
electric current generated. However, the chiral charge difference will also get contributions 
due to non-trivial vorticity and background magnetic field. Since, j M = nu M + + ^bB^, 
the contribution due to both vorticity and magnetic field is 



This chiral charge difference will induce an electric field resulting in a combination of chiral 
magnetic effect and chiral vortical effect. In the case of plasma in early universe, the above 
expression will be proportional to the rate of baryogenesis. 

6 Conclusion and future directions 

Chiral magnetic and vortical effect due to sphalerons has evinced much interest recently 
because it is a candidate model for explaining the observed charge dependent azimuthal 
asymmetries in the heavy ion collisions. [31] We attempted to model the charge asymmetry 
generation process in quark gluon plasma by constructing in this work explicit, analytic 
solutions of relativistic hydrodynamics containing parity violating and anomalous terms. 
Some solutions constructed in sections 3 and 4 are devoid of electromagnetic fields, though 
they possess non-trivial vorticity. These can nevertheless be candidate model for processes 
involving chiral vortical effect. The sphaleron solution constructed in section 5 is richer and 
have non-trivial values for all the parity violating and anomalous terms possible at the first 
order in the hydrodynamic approximation. It displays a combination of both chiral vortical 
and chiral magnetic effect. Azimuthal correlations between 3 particles are used to measure 
charge separation effects in experiments. They have been shown to be related to correlators 
of coefficients of sine terms in the Fourier decomposition of the charged particle azimuthal 
distribution. [32] Such sine terms are manifestly parity violating. The conjectured forma- 
tion of metastable topological configurations is in line with observed differences between 
same charge and opposite charge correlators. These effects disappear at low energies which 
is consistent with the sphaleron model. [33] The equations of hydrodynamics admit many 
different solutions and it is a priori difficult to say whether a particular solution actually 
occurs or not. However, the charge separation can be due to many such sphalerons occur- 
ring in different point of space and time in the quark gluon plasma. Further observations 



An 




(5.12) 



Hence, the total chiral charge difference created is 




(5.13) 
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are needed to fully ascertain whether these parity violating effects can be attributed to such 
configurations. In this direction, a detailed study of such solutions may facilitate in figuring 
out some typical predictions out of them and in providing theoretical predictions/hints for 
future tests. 

We believe our solutions also to be relevant in many other different contexts like plasma 
in early universe, superfluids and Fermi liquids as discussed in the introduction. Another 
interesting application area of our solutions can be neutron stars. The core of a neutron star 
is made up of highly dense quark matter displaying superfluid properties. [34, 35, 36] We 
anticipate that our results will also be relevant in describing some properties of vortices in 
neutron star. It will be interesting to construct explicit models based on these solutions to 
quantitatively demonstrate and evaluate the significance of the chiral magnetic and chiral 
vortical effect in various situations. It will help us to chart out the kinetics and dynamics 
of many processes in greater detail. 

We will also like to see the role played by dissipative coefficients in these solutions 
and in the processes in general. Quark gluon plasma possesses negligible viscosity and 
it is considered nearly perfect fluid. However, inclusion of dissipation in these solutions 
will enable us to make more realistic models and make better quantitative predictions out 
of them. Finally, from fluid dynamic point of view, the set of equations given by eqs. 
(4.10) is one of the main outcomes of our analysis. These equations can be used to gen- 
erate a variety of solutions of relativistic hydrodynamics relevant in different contexts. It 
is also an interesting problem to investigate that what kind of fluids i.e. fluids specified 
by equation of state can accommodate hydrodynamic solutions governed by the reduced 
set of equations written in eq. (4.10). We will also like to check the stability and linear 
response of our solutions against various perturbations. Also, we have only demonstrated 
that steady state solutions of eq. (4.10) are available, but it will be interesting to construct 
their explicit time dependent solutions or to demonstrate the consistency of these equa- 
tions including the temporal derivative terms. Nevertheless, our solutions being explicit 
analytic solutions, can also be used as benchmarks to check numerical codes for relativistic 
hydrodynamics. Analytic solutions of non-relativistic hydrodynamics have long been used 
for such purposes. 
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